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Abstract. There has been continued numerical work on the identification of scaling 
dimensions of an integrable two-dimensional lattice model of polymers via transfer 
matrix methods. The polymer model is known as the semi-flexible VISAW model 
and the integrable point is located at a particular set of Boltzmann weights. This 
point has been argued to describe the collapse phase transition of the polymer. We 
directly investigate the polymer scaling exponents, related to the scaling dimensions, 
via Monte Carlo simulations using the PERM algorithm. The polymer model was 
originally investigated via the exact solution of the 0{n) loop model and the polymer 
exponents, v = 12/23 ~ 0.522 and 7 = 53/46 ~ 1.152, were subsequently found via 
identification of the scaling dimensions x t = 1/12 and Xh = —5/48. By simulating this 
polymer model for walks up to length 4096 we find v = 0.576(6) and 7 = 1.045(5), 
which are clearly different from the predicted values. Our estimate for the exponent v 
is compatible with the known 0-point value of 4/7 and in agreement with very recent 
numerical evaluation by Foster and Pinettes pQ 
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1. Introduction 

The study of the critical properties of lattice polymers, and thus of 0(n) models when 
we let n — > 0, in two dimensions has been ongoing over decades theoretically and 
numerically. Nienhuis in 1982 [2] was able to compute the critical exponents of free self- 
avoiding walks (SAW) by considering a model of non-intersecting loops on the hexagonal 
lattice, and in 1987 Duplantier and Saleur [3] were able to model bond interactions 
introducing vacancies on the same lattice, obtaining a full set of critical exponents for 
the polymer collapse transition, a tri-critical point, in the interacting self- avoiding walk 
(ISAW) model. Their conjectured values for the exponents have been subsequently 
confirmed numerically by Prellberg and Owczarek in 1994 |4J. We will refer to (the 
universality class of) this critical point as the '#-point'. 

In the quest for a solvable 0(n) model on the square lattice, Blote and Nienhuis in 
1989 [5] considered a lattice model (related to the Izergin-Korepin vertex model), which 
includes energies for site-collisions and straight segments. For this model five critical 
branches are exactly known [5j El [7J |8] . One of these branches (named 'branch 0' in [5]) 
is similar to the #-point previously found by Duplantier and Saleur [3] (for this branch 
Batchelor j9] obtained the exponents i/ = 4/7, 7 = 6/7), and two branches correspond 
to dense and dilute (SAW) polymers. The two remaining branches are, respectively, 
associated with a combination of Ising-like and 0(n) critical behaviour and with a new 
tri-critical point. This other tri-critical point (with exponents v = 12/23, 7 = 53/46) 
[S] is another candidate for describing a collapsing polymer. The polymer model [TU] 
associated with this particular 0(n) model is a variant of the vertex interacting self- 
avoiding walk model (VISAW): one adds a stiffness parameter and at the same time 
excludes so-called crossings. The resulting set of configurations is actually a subset 
of self-avoiding trails (bond-avoiding walks). The integrable point found by Blote and 
Nienhuis, which we shall refer as a the BN-point, corresponds to a particular value of 
the stiffness and interaction weights. When stiffness vanishes [llj #-point behaviour 
ensues. 

Foster and Pinettes {TOj have studied the semi-flexible VISAW at this special BN- 
point, and have also studied the VISAW model without stiffness, using the corner 
transfer matrix renomalisation group method. Some agreement and some discrepancy 
with the scaling dimensions proposed [3 El [8] was found in [10] and a first order 
nature to the transition was conjectured. Very recently Foster and Pinettes pQ have used 
transfer matrices and the Density Matrix Renormalisation Group Method (DMRG) to 
consider the bulk and surface exponents of these models. They have found values of the 
exponent v much closer to 4/7 than 12/23. In this paper we study by means of Monte 
Carlo simulation the semi-flexible VISAW polymer model precisely at the BN-point. We 
find estimates for the exponents, and hence the scaling dimensions, that are in harmony 
with those found by Foster and Pinettes |T] and at variance with those predicted by 
Warnaar et al. [S] . 
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2. Semi-flexible VISAW 



The semi-flexible VISAW model can be denned as follows. Self-avoiding trails, or simply 
trails, are lattice paths that can be formed such that they never visit the same bond 
more than once. Such paths can generally visit the same site of the lattice either by a 
collision, where the trail touches itself, or via a crossing, where two straight segments 
of the path cross over one another. Consider the subset of bond-avoiding lattice paths 
(trails) on the square lattice, V n , where no crossings are allowed. Given such a restricted 
trail ip n € V n , we associate an energy — e t every time the path visits the same site more 
than once, which it can only do by colliding with itself: see Figure [TJ Additionally, we 
define a straight segment of the trail by two consecutive parallel edges, and we associate 
an energy — e s to each straight segment of the trail. 
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Figure 1. An example of a semi-flexible VISAW configuration with three (m = 3) 
collisions associated with a Boltzmann weight r and five (s = 5) straight segments, 
each associated with a Boltzmann weight p. 



For each configuration ip n e V n we count the number m(ip n ) of doubly-visited sites 
and s(ip n ) of straight segments: see Figure [TJ Hence we associate with each configuration 
a Boltzmann weight T m (i>n)ps(ip„) w h ere T = exp(/3e t ), p = exp((3e s ), and (3 is the inverse 
temperature l/fc^T. The partition function of the model is given by 

Z n (r,p) = T m ^ n) p s ^ n \ (2.1) 

V>neV n 

The finite-length reduced free energy is 

k„(T) = -log Z n (2.2) 
n 

and the thermodynamic limit is obtained by taking the limit of large n, i.e., 

k(T) = lim K n (T). (2.3) 

n— >oo 

It is expected that there is a collapse phase transition at a temperature T c characterized 
by a non-analyticity in k(T). Equivalently, one can think of varying r at fixed p so that 
there is a collapse at some value of r = r c (p). 

The probability of a configuration ip n is then 

p{^r,p)= f- , (2.4) 
Zn{r,p) 
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and the average of any quantity Q over the ensemble set of path V„ is given generically 
by 

(Q)n(r,p)= Q^n)p^n,r,p). (2.5) 

In this paper we are interested in the following quantities. We calculate three 
measures of the size of the polymer, (R%) n , (Rm)n and (R^) n , defined as follows. We 
specify any n-step path ip n on a lattice by a sequence r , rj., . . . , r n of vector positions 
of the vertices of that path. Firstly, we are interested in the average-square end-to-end 
distance 

(Rl) n = (r n • r n > , (2.6) 

secondly, the ensemble average of the mean-square distance of a monomer from the 
endpoints 

1 - 

(^>»=- T £< r *- r <>' (2-7) 
and defining the average centre-of-mass as 

1 n n 

V ' i=0 j=0 

we, thirdly, are interested in the average radius-of-gyration 

(R 2 9 )n = (R 2 Jn - (RX ■ (2.9) 

In the above formulae we use r = 0. 



2.1. Scaling 

The partition function at, and above, the collapse temperature is believed to scale as 

Z n ~ DfMV 1 (2.10) 

where \x is known as the connective 'constant' and is related to the thermodynamic free 
energy via 

/x = e K(T) . (2.11) 

The constant D is also temperature dependent but 7 is expected to be universal, 
depending only on the temperature in as much as its value is above, or at, the collapse 
temperature. In two dimensions it is well established [2] that for T > T c (equivalently 
r < T c (p)) we have 7 = 43/32. 

The collapse transition can also be characterized via a change in the scaling of the 
size of the polymer with temperature. The three measures of the size of the polymer 
defined above are expected to scale as 



Rl ~ C R n 2 > 



(2.12) 
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where the amplitude Cr is non-universal and temperature dependent, while v is expected 
to be universal, depending only on the temperature in as much as its value is above, at, 
or below, the collapse point. In two dimensions it is also estabslished j2] that v = 3/4 
for T > T c . 

Duplantier and Saleur [3] identified a tri-critical point, known as the #-point, which 
is expected to describe the collapse of a polymer in two dimensions. This point has 
thermal and magnetic scaling dimensions x t = 1/4, Xh = and consequently polymer 
exponents v = 4/7 and 7 = 8/7. On the other hand, Warnaar et al. [8] predicted 
v = 12/23 and 7 = 53/46 for the semi-flexible VIS AW model at its collapse point. 



2.2. Amplitudes 

One can also usefully define the finite-length amplitude ratios 



An = 7xt and B n = { -^, (2.13) 



(RDn m 

since these approach universal values [12 



A n ^A OQ = -^- and B n ^B 00 = < ^ (2.14) 

in the limit n — > 00. For collapsing polymers, the limiting values should depend only 
on dimension and whether the temperature is above or at the collapse transition point. 

For free self-avoiding walks (which should include the VISAW model at high 
temperatures) it was predicted [131 EE] that 

XA^ - IB^ + - = . (2.15) 

In the derivation [13] of this invariant the factor multiplying was given by 

A = 2 + —, (2.16) 

Vh 

where y t = 4/3 and = 91/48 are the thermal and magnetic renormalisation group 
eigenvalues, respectively, of the dilute O(0) model. These eigenvalues are related to the 
conformal scaling dimensions via y = 2 — x. Hence 

X( XuXh ) = 2+l—^ . (2.17) 
2 -x h 



The identity (2.15) implies that one can estimate this function of the scaling dimensions 

as 

4B — 1 

X(x t ,x h )= ^ , (2.18) 

from estimates of A^ and B^. This was done for various collapse models in [15j. Hence, 
if we have a conjectured value of one of the scaling dimensions we can estimate the other 
from an estimate of A. 
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3. Integrable Blote-Nienhius Point 

The special multi-critical point of the 0(n) model that maps to the semi-flexible VIS AW 
and allows for the calculation of the scaling dimensions via the Bethe Ansatz is given 
by special values of the parameters in the grand canonical partition function 

oo 

G(K;r, P ) = J2K n Z n (r,p). (3.1) 

n=0 

The location of this point is reported exactly in [5] 

w = Klr bn = {2 - [1 - 2sin(0/2)][l + 2 sin^/2)] 2 }" 1 
K bn = -Aw sin(0/2) cos(tt/4 - 0/4) 
p bn K bn = w{l + 2sin(6/2)} 

9 = —7i /A (branch 3 in [5]) . 

Alternatively, this can be expressed in explicit algebraic numbers or evaluated 
numerically as K bn = 0.446933 ...,p bn = 0.275899 . . ., and r bn = 2.630986 .... 

From the value of K bn we can give the conjectured value of the connective constant 



\i in the canonical model (2.1) at the Boltzmann weights (r,p) = (tbn,Pbn) as 



Knn,Pbn) = ^-= ^ + ^V 2 "^2+ \/ 2 -^J V2-V2-V2 

= 2.23746994... (3.2) 
The set of scaling dimensions evaluated in [H] at this (multi-)critical point are 

xi = ^ - \ for leN. (3.3) 
16 6 

The "thermal" x t and "magnetic" Xh scaling dimensions were identified [S] as 

5 1 

Xh = xi = -— = -0.10416. . . and x t = x 2 = — = 0.08333 ... . (3.4) 
48 12 

It can be seen that these scaling dimensions are not those of unweighted self-avoiding 

walks; this is compatible with the hypothesis that they are those of a collapse multi- 
critical point. 

The exponents v and 7 were then calculated [8] in the standard way as 

1 12 . , 

0.52174... (3.5) 



2 - x t 23 

53 

7 = 2v (1 - x h ) = — = 1.15217 .... (3.6) 
46 
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4. Simulation results 



We simulated the model (2.1) at the fixed values of p = pb n and r = r bn using the 
PERM algorithm. We ran three simulations with maximal length N max = 1024, 2048 
and 4096, growing S ~ 10 7 independent walks each and collecting from 8.7- 10 7 to 1.5-10 8 
samples at each maximal length. The number of samples adjusted by the number of 
their independent growth steps is between 2.1 • 10 5 and 4.8 • 10 5 "effective samples". 




Estimates of v 
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Figure 2. A double-logarithmic plot of three different measures R\ of the polymer size, 
the mean-square end-to-end distance {R%)n (top), the square of the radius of gyration 
(Rg) n (bottom), and the mean-square distance of a monomer to its end points (R^ n ) n 
(middle) versus the length n of the path. The data from different simulations up to 
lengths 1024 (blue circles), 2048 (green triangles), and 4096 (red line) clearly overlap. 
The exponent estimates for v shown are derived from fitting straight lines through 
lengths 1024 to 2048 and 2048 to 4096, respectively. 



In Figure [2] we plot on a double-logarithmic scale the three different measures R 2 n of 
the polymer size. From various fits we consistently find estimates of v near 4/7 rather 
than 12/23. Our best estimate is 



v = 0.576(6) . 

This leads to an estimate of the thermal scaling dimension as 



(4.1) 



x t = 0.26(2) . 



(4.2) 



To obtain an estimate of 7 we looked at the scaling of the canonical partition 
function (2.10). We first measured fi by a simple linear fit obtaining 2.2375(1), and 
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Figure 3. A double-logarithmic plot of the normalized partition function Z n / fi n versus 
the length n of the path. The difference between data from different simulations up to 
lengths 1024, 2048, and 4096 indicate the accuracy achieved. The exponent estimate 
for 7 shown is derived from fitting straight lines through lengths 256 to 1024 to data 
from the best-converged simulation. 



observing that our value matches the value obtained from the BN-model, we then 
assumed /i = /i& n . We hence plotted, in Figure [3j on a double-logarithmic scale the 
normalized partition function Z n /n^ n versus the length n of the path. 
This allows us to estimate 7 from straight line fits which we give as 

7 = 1.045(5) . (4.3) 

Interestingly, this value is different from both the #-point value of 8/7 = 1.4228 . . . 
and the BN-point 53/46 = 1.1521 . . .. Using 7 = 2z/(l — Xh) and our estimates of 7 and 



v in equations (4.3) and (4.1) gives us the estimate 

x h = 0.093(13) . (4.4) 
We point out that this estimate is positive while the conjectured value above in 



equation (3.4) is not. 

To obtain an independent estimate of Xh we attempted to estimate the universal 
quantity A described above. In Figure [4] we plot the finite-size amplitude ratio 
combination A n = (4S n — l)/(2A n ) versus n~ 4 / 7 , which is the natural scale given the 
results above for the size measures. We find an estimate of this universal value as 



A = 2.93(3) 



(4.5) 
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Figure 4. A plot of the finite-size amplitude ratio combination A n versus n~ 4 / 7 , 
together with a straight-line fit. The difference between data from different simulations 
up to lengths 1024, 2048, and 4096 indicate the accuracy achieved. The extrapolated 
value for A shown is derived from fitting straight lines through lengths 256 to 1024 to 
data from the best-converged simulation. 



Unfortunately, the error estimate here is relatively larger than that estimated from the 
partition function analysis and we estimated Xh = 0.12(4), which encompasses our more 



precise estimate in (4.4). 

For the sake of completeness we have also found estimates of the universal amplitude 
ratios A and B by fitting against a correction of n~ 4 ^ , which provide consistent straight 
fit extrapolations 

A = 0.1534(10) and B = 0.475(5) . (4.6) 

These values of A and B are different from the values for any of the three collapse 
models considered in 1151. 



5. Conclusions 

We have simulated the special point, known as the BN (Blote-Nienhuis) point, of the 
semi-flexible VISAW model of polymer collapse, which is associated with an integrable 
branch of the 0(n) loop model [5]. The exponent estimates we find, v = 0.576(6) 
and 7 = 1.045(5), are not in accord with those previously found from the 0(ri) loop 
model. Our estimate of v is compatible with the value accepted for the ^-point, which is 
4/7 = 0.5714 . . ., and in good agreement with the results of Foster and Pinettes pQ who 
have used transfer matrices and DMRG. However, our estimate of 7 is not comparable to 
any known value. We have found estimates of 7 via two different methods. One method 
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we used involved the direct estimation of the exponent from the partition function, and 
the other used results from conformal field theory and universal amplitude ratios of 
different size measures of the polymer; our estimates from these two methods broadly 
agree. Clearly, something subtle is occurring in this system if our numerical analysis 
is accurate. Of course, large corrections to scaling may be at work here. In any case, 
further theoretical work is needed to tease out this issue. 
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